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Abstract 

We study how the determination of the Hubble constant from cosmological dis- 
tance measures is affected by models of dark energy and vice versa. For this purpose, 
constraints on the Hubble constant and dark energy are investigated using the cos- 
mological observations of cosmic microwave background, baryon acoustic oscillations 
and type la suprenovae. When one investigates dark energy, the Hubble constant 
is often a nuisance parameter, thus it is usually marginalized over. On the other 
hand, when one focuses on the Hubble constant, simple dark energy models such as 
a cosmological constant and a constant equation of state are usually assumed. Since 
we do not know the nature of dark energy yet, it is interesting to investigate the 
Hubble constant assuming some types of dark energy and see to what extent the 
constraint on the Hubble constant is affected by the assumption concerning dark 
energy. We show that the constraint on the Hubble constant is not affected much 
by the assumption for dark energy. We furthermore show that this holds true even 
if we remove the assumption that the universe is flat. We also discuss how the prior 
on the Hubble constant affects the constraints on dark energy and/or the curvature 
of the universe. 



1 Introduction 



Recent precise cosmological observations provide us with a large amount of data to probe 
the evolution and the present state of the universe. We now customarily extract the 
information from them by constraining cosmological parameters. The Hubble constant, 
the expansion rate of the universe at present, is one of the most important cosmological 
parameters and measured in various ways. 

In most studies which consider the constraint on the Hubble constant from cosmological 
observations, a cosmological constant is often assumed as dark energy. However, the 
nature of dark energy is not understood yet, which is one of the challenging problem in 
cosmology today, thus the constraints on the Hubble constant should also be considered 
assuming some possible dark energy models other than a cosmological constant. Since the 
energy density of dark energy can change some distance measures which have been used to 
constrain the Hubble constant, the assumption concerning dark energy is expected to have 
an influence on the constraint on it. In some works, the constraint on the Hubble constant 
is obtained without assuming a cosmological constant (e.g. Refs. [1,2]) but a constant 
equation of state is still assumed. However, a lot of models of dark energy proposed 
so far have a time-varying equation of state and most of recent works of dark energy 
accommodate its time dependence in some way. Thus, in light of these considerations and 
the importance of the Hubble constant, it is interesting to study constraints on the Hubble 
constant including timc-cvolving dark energy equations of state. This is one of the issues 
which we are going to investigate in this paper. 

On the other hand, notice that the Hubble constant is often treated as a nuisance 
parameter over which we marginalize when we investigate dark energy since the parameters 
for dark energy themselves are those of interest in such a case. However, it should be 
mentioned that, when one studies the nature of dark energy, observations of distance 
measures such as the angular diameter distance which is relevant to the position of acoustic 
peaks in cosmic microwave background (CMB) power spectrum and the scale of baryon 
acoustic oscillations (BAO) are often used. Since the Hubble constant affects such distance 
measures, it is expected that the prior on the Hubble constant can affect constraints on 
the nature of dark energy. Furthermore, the Hubble constant can also be determined with 
the cosmic distance ladder measurements independently from the cosmological distance 
measurements mentioned above. Hence investigating dark energy with some priors on the 
Hubble constant is also an interesting subject, which is discussed in this paper too. 

To study the issues mentioned above, we study the constraints on the Hubble constant 
and dark energy using the observations of CMB, BAO and type la supernovae (SN). For 
dark energy, we assume some types of time-varying equation of state as well as the case 
with a constant equation of state. By assuming several priors on the Hubble constant and 
dark energy equation of state, we can see how the determination of one of them can affect 
that of the other. In addition, we study them allowing a non-flat universe to make our 
analysis more general. 

The structure of this paper is as follows. In the next section, we summarize our method 
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for obtaining constraints from cosmological observations. Some parametrizations of dark 
energy adopted in this paper are also briefly explained. The data used for the analysis 
are mentioned there too. In Sec. [3} we present our results and discuss some implications 
of our results for the study of the Hubble constant, dark energy and the curvature of the 
universe. In the final section, we summarize our results and give the conclusion. In the 
appendix, we give some quantitative explanations how the combinations of CMB, BAO 
and SN can break degeneracies among the Hubble constant, matter density, dark energy 
parameters and the curvature of the universe. 



2 Dark Energy Parametrizations and Method 

In this section, we explain the method to constrain the cosmological parameters such as 
the Hubble constant and dark energy equation of state. We use the data from CMB, 
BAO and SN. Before we describe the method in detail, first we mention the dark energy 
parametrizations adopted in this paper, which accommodate the time variations of its 
equation of state. Among various possible parametrizations of dark energy equation of 
state, we use simple and often used ones. The first parametrization adopted in the following 
analysis is [3,4]: 

2^ 

Wx=WQ^{^-a)W\ = WQ^ ——w\, (1) 

which varies in proportion to the scale factor a (normalized as a = 1 at the present epoch). 
If we define dark energy as a fluid which can accelerate the universe, this can motivate 
the following prior: 

< (2) 
For the particular parametrization Eq. ([T|, this is ensured by forcing wq and wi to satisfy 

11 

wo + wi<--, wo<--. (3) 

In principle, a dark energy fluid can have its equation of state larger than —1/3. Thus a 
weaker prior wx < is sometimes adopted, requiring that dark energy does not dominate 
the universe at early time, say at the epoch of recombination, since such early domination 
of dark energy is apparently inconsistent with cosmological data. However this prior allows 
the case with wx ~ where dark energy behaves in almost the same way as dark matter as 
far as the background evolution is concerned. In such a case, dark energy can be considered 
as a part of dark matter, which complicates the interpretation of the constraints on dark 
With the prior Eq. (|2|, the dark energy cannot be a dominant component of 
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This holds true when one considers the background evolution alone, which is the case of the present 
paper. If we include the information of perturbation of dark energy, the effect of dark energy can be 
distinguishable even if its equation of state is almost the same as that of dark matter. In that case, we 
need a full CMB angular power spectrum analysis to obtain the constraint from CMB but such analysis 
is beyond the scope of this paper. 
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the universe at early times and even its fraction to the total energy density is negligible. 
Thus we can safely avoid the complication mentioned above by using the prior Eq. g. We 
provide more detailed discussion on the effects of dark energy perturbation on the CMB 
power spectrum in Appendix [A| 

We also adopt another type of dark energy parametrization: 

Wo H z (for z < z^) 

(4) 

(for z > z^:), 

where we interpolate the value of wx linearly with respect to the redshift z from the present 
epoch back to some transition redshift z^,. wx becomes wq at z = and wi for z > z^. 
The case with wq = wi = —1 corresponds to a cosmo logical constant. The prior Eq. ^ 
is satisfied by setting wq < —1/3 and Wi < —1/3. This parametrization is essentially the 
same as the one which has a linear dependence on the redshift z such as wx = Wa+wi, z. In 
fact, this form is adopted in many literatures with a cutoff at some redshift to avoid large 
value of Wx at z ^ 1000 which is relevant to the CMB constraint [5-7]. In the following, 
we consider several values of z^, to see how the dark energy parametrization affects the 
constraints on the Hubble constant and vice versa. We note that this parametrization 
assumes a finite derivative of wx at the present epoch so we cannot probe the dark energy 
model which behaves like a cosmological constant (or a constant wx) from some earlier 
epoch to the present but variable before that epoch. Also, we should note that although 
this parametrization appears to have three parameters, we regard it as a family of two- 
parameter models labeled by z^:. Namely, we consider several two-parameter models each 
with a different value of z*. Such a distinction among dark energy parameters may be 
artificial, but it is simpler and sufficient for our goal to show how priors for the Hubble 
constant and ways of parametrizing the dark energy equation of state (and the assumption 
for the curvature of the universe) could influence a conclusion drawn from the analysis of 
CMB, BAG and SN data. 

Now we give detailed descriptions of the method how we make use of observational data. 
Incidentally, the model parameter fitting to the data is performed by minimization. We 
use the Brent method [8] extended to multi-parameters, as described in Ref. [9], to search 
a minimum efficiently. 

2.1 CMB 

To fit a model to the CMB data, we first note that we here do not use the whole information 
of the CMB power spectrum but use only the acoustic scale 9a and the matter density Um 



as Eq. ( 18 ) *^ , This is because if we include the effects of perturbation of dark energy and 



'^^ Recently, Refs. [10, 11] noted that it is better to use = tt/6a and R = ^J^l^^Jr[^rg{z-ccc) simultane- 
ously to obtain the constraint from CMB. We in this paper use rather conventional way of using 9a and 
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perform a full CMB angular power spectrum analysis to obtain the constraint from CMB, 
it would be very time-consuming. One of the purposes of this paper is to investigate 
cosmological constraints assuming various dark energy parametrizations. Thus, in this 
case, much faster method is preferable and it is known that if dark energy dominates the 
universe only at late time, the constraint derived from the information on the background 
evolution well captures the nature of dark energy. This condition is satisfied by adopting 
the prior Eq. (|2|. 

The acoustic scale 6a which defines the characteristic angular scale of the acoustic 
oscillations is written as 

Oa = (5) 

9a is given once we determine the comoving angular diameter distance to the last scattering 
surface re{zrec) and the sound horizon at the recombination epoch rs(zrec) where z^ec is the 
redshift of the epoch of recombination. The comoving angular diameter distance to the 
last scattering surface is given as 

^ jj^^ {^^r w) ■ 

where S is defined as S{x) = sin(x) for a closed universe, S{x) = sinh(x) for an open 
universe and S{x) = x with the factor A/|fifc| being removed for a fiat universe. Hq is the 
Hubble constant and we sometimes denote it using a renormalized quantity h defined as 
Hq = lOO/i km s~^Mpc~^. The expansion rate H{z) is given as 



H\z) = Hi 



fi,(l + ;z)^ + fi„(l + ;z)' + fifc(l + z)" + fixexp ( 3 / (l + Wx{z)) 







(7) 

where f2r,fim,fifc and fix represent the present values of the energy density of radiation, 
matter, curvature and dark energy normalized by the critical density. Note that fi^ = 
1 — fir — fim — fix- The sound horizon at recombination is 

f^'"" Cs 

rs{z^ec) = / -T¥?da (8) 



/o am 

where a^ec = 1/(1 + ^rec) and Cg is the sound speed of the photon-baryon fluid: 

Pj + Pb 3(1 + 



ujm as in, for example, Refs. [12,13], but both ways make use of the same features of the CMB power 
spectrum, the peak position and height, and should give similar result. We checked that the constraints 
on the ^lm~h plane from these two methods axe almost same for the case with a cosmological constant. 
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with R — Spb/Ap^ being the scale factor normahzed to 3/4 at the photon-baryon equahty. 
Since we consider the case where dark energy is irrelevant at the epoch of recombination, 
the above expression can be integrated analytically as [14, 15] 

rs(^rec) = -===ln<^ \ Mpc, (10) 

where Ub and Um are normalized baryon density and matter density (i.e., Ui = VLih?). The 
radiation-to-matter and baryon-to-photon ratios at the epoch of recombination. Tree and 
it!rec) are given as 

r,ec = 0.042 o;-^(WlO'), (H) 
i?,ec = 30o;b(WlO=')"'- (12) 

In fact, the redshift at the epoch of recombination slightly depends on the energy density 
of baryon and matter. We include the dependence by adopting the fitting formula z^-^^ [16]: 

Zrec = 1048 [1 + 0.00124 [1 + <] , (13) 

where the functions gi and g2 are given as 

g^ = 0.0783 [i + 39.5 u;^-^^^] -i ^ (^4) 

^2 = 0.560 [1 + 21.1 (15) 

We fix the value of uj^ with the mean value from the WMAP3 analysis uj^ = 0.02229 which 
is obtained for the ACDM model. 

Once we give the Hubble constant, matter density, baryon density and dark energy 
parameters, 9a can be calculated. To constrain the cosmological parameters, we use the 
value of Oa reported by WMAP3 for ACDM 

OA,ohs = 0.5952 ° ± 0.0021 °. (16) 

In addition to 9a, we also use the prior for ujm which is given by WMAP3 for ACDM as 

c^m,obs = 0.1277 ±0.008. (17) 

Notice that the values of ^A,obs and a;^^obs are almost unchanged even if we consider a 
model with a constant equation of state or a non-flat universe [1, 11]. 
Using 9a and o;^, we calculate the from CMB as 

2 (^A — ^A,obs)^ , (<^m — <^m,obs)^ o\ 
XCMB - ^ \ ^ ■ lloj 

Roughly speaking, in this way, although we are not using the whole information of the 
CMB power spectrum but using only the information of the first peak position and height 
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(after we extracted Ub from the second peak height to calculate 9a), since 9a and Um are 
determined from such different, horizontal and vertical, features of the CMB power spectra, 
the fo^^ the parameters we concern in this paper can be constructed as above. This is 
explicitly justified in Sec. |3]for the fiat ACDM model. Furthermore, since the effects of the 
curvature of the universe and dark energy equation of state on the CMB power spectrum 
mainly appear as the change of 9a , we can use this Xcmb constraints involving these 
parameters. 



2.2 BAO 

We also use the baryon acoustic oscillation scale measured by SDSS [25]. To take the BAO 
data into account, we use the parameter Dy which is defined as [25] 



Dviz) 



Hiz) 



1/3 



(19) 



where rQ{z) is the comoving angular diameter distance defined in the previous subsection*^ 
This quantity depends on the scalar spectral index and oji, slightly. We take Ug = 0.958 
and Ub = 0.02229 which are the mean values for the ACDM model from WMAP3 data 
alone and calculate 1^^(0.35) following the procedure of Ref. [25] to find Dv{0.35) = 
1402 Mpc. Then we use 

^^(0.35)0^ = 1402 ± 64 Mpc, (20) 
for the constraint from BAO. The value of is calculated as 

2 _ {Dv{0.35) - Dvi0.35U,f 

xbao — ~2 ■ y-^^) 

^Dv{0.35) 

We note that the measurement of (0.35) assumes a flat ACDM model and there is no 
explicit check that the value does not change in a non-flat model and/or dark energy models 
with wx 7^—1- However, since the effect of non- flatness or dark energy is considered to 
be only the geometrical one (their effects on the perturbation do not affect the peak of 



'^■^ The curvature of the universe and dark energy equation of state also affect CMB power spectrum 
on large scales through the late-time integrated Sachs- Wolfe effect, but since the error from the cosmic 
variance dominates on large scales, such an effect can be neglected in most cases [17,18] . On the other hand, 
it should be mentioned that such large scale fluctuation may be interesting when one considers scenarios 
such as dark energy isocurvature fluctuation [f9, 20], the nature of dark energy perturbation [f8,2f-24] 
and so on. 

'^'^ In most of works of dark energy, the so-called A parameter is used. However, we use Dy{z) which 
explicitly depends on the Hubble constant for the purpose of this paper. The relation between Dv{z) and 
the A parameter is 

A X 0.35 „ „ /A X 0.35' 



DviQ.'i^) = ' ^ = 3.0 X 10-^ ( — ) Mpc. 
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the galaxy correlation function which is located at smaller scales), we can use Dv{0.35) 
to probe the background evolution just as we use 9a of CMB if the cosmology is not 
radically different from a flat ACDM model. This condition is considered to be satisfied 
in our analysis because we place the prior Eq. ^ discussed in the previous section and 
our resulting constraints as shown later are not far away from a fiat ACDM model when 
we combine with CMB and SN. 

2.3 SN 

As for SN data, we calculate the distance modulus 

/i = m-M = 51ogrfL + 25, (22) 

where m is apparent magnitude and M is absolute magnitude. Here di is the luminosity 
distance in units of Mpc which is written as 




Although the luminosity distance depends on the Hubble constant explicitly, its depen- 
dence is indistinguishable with the uncertainty in M. Thus when we calculate x^? the 
dependence on h vanishes by marginalizing over M as a nuisance parameter. 

For the analysis, we use 182 SNe from the Gold data set of Riess et al. (Gold06) [26] 
or 192 SNe of Davis et al. (Davis07) [26-28]. In this paper, we present our results for 
these two data sets separately. (Namely we do not combine these two data sets.) 

The main difference between these two data sets is that Davis07 uses the first data 
release of the ESSENCE supernova survey [27,29] at the intermediate redshift. In detail, 
Gold06 consists of 119 SNe (38 SNe are nearby, z^O.05) from the previous Gold data 
set [30], 16 SNe which are recently discovered by Hubble Space Telescope (HST) [26] and 
47 SNe from the first data release of the SNLS project [31]. Davis07 includes 60 SNe 
from ESSENCE, 57 SNe from SNLS, 16 SNe of Ref. [26], 14 SNe discovered by HST 
which are included in the previous Gold data set and 45 nearby (^^0.05) SNe. Thus, 
30 SNe which are discovered by HST and have relatively high redshifts are common in 
both data sets. Also there are 25 nearby SNe which they have in common. The SNLS 
first data release [31] has 73 SNe but Gold06 and Davis07 place a different criterion for 
the selection. We may roughly say that Davis07 is constructed from Gold06 by replacing 
SNe data at the intermediate redshifts which are discovered in earlier times by the recent 
ESSENCE data. 

2.4 Hubble constant 

One of the purposes of this paper is to constrain the Hubble constant using the cosmological 
observations introduced above. However, as mentioned in the introduction, the Hubble 
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Refs. 


Ho ± (statistic) ± (systematic) 


Freedman et al. [32] 
Sandage et al. [33] 
Macri et al. [34] 


72 ± 3 ± 7 
62.3 ± 1.3 ±5.0 
74 ± 3 ± 6 



Table 1: The values of the Hubble constant in unit of kms ^Mpc ^ determined from 
cosmic distance ladder measurements. 



constant can be measured with other methods such as using the cosmic distance ladder 
measurements. We summarize the recent values obtained with the distance ladder in 
Table [11 For a recent review on the Hubble constant, see e.g. Ref. [35] We would like 



to compare such Hq values measured in a relatively direct way with those obtained from 
the measurements of CMB, BAO and SN. Furthermore we also investigate constraints on 
dark energy putting a prior on the Hubble constant to see how the determination of the 
Hubble parameter affects the constraints on dark energy. Thus we also use these values as 
priors when deriving constraints on dark energy parameters and/or the curvature of the 
universe. When we use the prior on the Hubble constant, we include them by calculating 
the which is given by 

Xk = ""'"f"""'' - (24) 

For the error ano, since we do not obtain a significant effect with those of the current mea- 
surements given in Table [l| as will be briefly discussed in the following section (Sec. 3.2) 



we take a hypothetical value = 2kms~^Mpc~^. This value is motivated by the ex- 
pected accuracy of 1% [34] which could be obtained through measuring maser distance to 
a large number of galaxies in the Hubble flow by planned radio telescopes such as Square 
Kilometer Array (SKA) [38]. We adopt a slightly more conservative value than that. 



In particular, interesting constraints on the Hubble constant are obtained via gravitational lens time 
delays and the Sunyaev-Zel'dovich effect, although the uncertainties are comparable or slightly larger 
than those of the distance ladder measurements. See Refs. [36] and [37] for the recent measurements by 
these methods. Since these two methods in principle depend also on a dark energy model, their future 
improvements are expected to be useful for probing the cosmological constraints on the Hubble constant 
and dark energy parameters independently from CMB, BAO and SN. 
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3 Results 



Now we discuss the constraints on the Hubble constant assuming some dark energy models 
including time-evolving dark energy equations of state described in the previous section. 
We also investigate how the constraint on dark energy parameters are affected by the 
assumption of the Hubble constant. Effects of the curvature of the universe in constraining 
these parameters and the constraint on itself are also investigated. 

3.1 Constraint on the Hubble constant 
3.1.1 Case with a cosmological constant 

First we show the constraints in the ^Im^h plane for the case with a cosmological constant 
m Fig.[lj Inthefi gure, we show the constraints from several data sets, i.e., the cases with 
CMS alone, CMB+BAO, CMB+SN and CMB+BAO+SN. We do not show the constraint 
from SN data alone and BAO data alone here. Since, for SN data, the dependence of 
the Hubble constant on the luminosity distance is absorbed into the uncertainty in the 
absolute magnitude of SN, which is marginalized over, we cannot obtain the information 
on Hq from SN alone. For the constraint from BAO, we use the parameter Dy for the 
analysis in this paper, which cannot constrain the Hubble constant much by itself from 
the currently available data. Thus we only show the constraints from CMB alone and the 
combinations of CMB and other data sets. 

For the analysis with CMB data alone, our method gives consistent results with those 
obtained by the WMAP team who uses information of the entire CMB power spectrum. 
For the marginalized la values of h and ilm, theirs are h = and flm = 0.241 ± 

0.034, which are satisfactorily close to our values h = 0.729±0.035 and Qm = 0.244±0.038. 
Closer comparison reveals that their two-dimensional contours in this plane (Fig. 10 in 
Ref. [1]) are somewhat larger than ours (Fig. [T] (a)) but the difference scarcely affects 
the estimation of h and Qrn as mentioned above. This small difference arises most likely 
because we fix Ub to the WMAP central value and our method cannot incorporate the 
uncertainty due to the value of Ug. As for the SN data sets, we show our results with 
Gold06 and Davis07 separately. In Table [2| we report the central values and la errors 
from various combinations of data sets for the case with a cosmological constant. The 
values of the minimum for each case are also shown. When we combine SN with 
CMB, the central value of h becomes smaller since the constraint on and h lies along 
the region around Qmh"^ being constant (for the case with WMAP3, flmh"^ = 0.1277) 
and SN data favor somewhat larger than CMB. This is more conspicuous for Gold06 
than Davis07. Another point is that the BAO data also favors a slightly smaller value of 
h. Hence when we use the data set of CMB+BAO+SN(Gold06), the favored value of h 
becomes considerably smaller. 
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Figure 1: Contours of la and 2a allowed regions in the ^mrh plane for the case with a 
cosmological constant in a flat universe, (a) CMB alone, (b) CMB+BAO, (c) CMB+SN 
and (d) CMB+BAO+SN. In the panels (c) and (d), we treat the SN data from Gold06 
(red solid line) and Davis07 (blue dashed line) separately, thus two different constraints 
are shown. 





h {la) 


Vtm {la) 


A mm 


^ of data 


CMB alone 


0.729 ± 0.035 


0.244 + 0.038 


0.00 


2 


CMB+BAO 


0.704 + 0.027 


0.270 + 0.033 


1.41 


3 


CMB+SN(Gold06) 


0.682 + 0.021 


0.299 + 0.028 


162.0 


184 


CMB+SN(Davis07) 


0.714 + 0.022 


0.259 + 0.025 


195.9 


194 


CMB+BAO+SN(Gold06) 


0.678 + 0.019 


0.304 + 0.026 


162.2 


185 


CMB+BAO+SN(Davis07) 


0.704 + 0.019 


0.269 + 0.023 


197.0 


195 



Table 2: The central values and la errors from several combinations of data sets for the 
case with a cosmological constant. The values of the minimum and the numbers of 
data are also shown. 
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3.1.2 Case with a constant equation of state 



Next we show the constraints for the case with a constant equation of state for dark energy 
(which corresponds to the cases with wi = in Eq. ([T]) and with wq = wi in Eq. Q). In 
Fig. |2| the constraints in the Qm^h plane are shown after marginalizing over wq. In the 
figure, 2a constraints from CMB alone (red solid line), CMB+BAO (green dashed line), 
CMB+SN (blue dotted line) and all combined (purple shaded region) are shown. We treat 
the SN data from Gold06 and Davis07 separately, thus we show the constraints using them 
in the left and right panels in Fig. |2] respectively. 

Since dark energy equation of state directly affects the angular diameter distance, the 
allowed region from CMB alone becomes significantly larger compared to the case with 
a cosmological constant, which means that there is a strong degeneracy among Qm, wq 
and h. However, if we combine other data set such as BAO and SN with CMB data, 
the allowed region becomes similar to that for the case with a cosmological constant (the 
region becomes larger but only a little). Although the preferred values of wq which is 
marginalized over in the figure are in general different for each observation for fixed h and 
Qm, only around the allowed region in the figure, those values from different observations 
happen to coincide to be Wq ~ —1 (i.e. a cosmological constant). Thus we obtain an almost 
same result as the cosmological constant case even for the case where we marginalize over 
a constant equation of state. We see that the combination of the data from CMB+SN 
seems to be enough to constrain h and flm whereas the CMB+BAO constraint is much 
weaker. This is because the degeneracy curves in the flm^wo plane extend to almost the 
same direction for CMB and BAO, but those of CMB and SN are complementary (see 
Fig. 6 in Ref. [39])*^ Thus as far as the degeneracy between flrn and Wq is concerned, the 
combination of CMB+SN gives a severe constraint. 



3.1.3 Case with time-evolving equations of state 

In Fig. [3} we show the results for the case with the time-evolving dark energy equation of 
state parametrized as Eq. ([T]). In the figure, the values of wq and wi are marginalized over. 
As the case with a constant equation of state, although the allowed region from CMB data 
alone is significantly larger compared to that for a cosmological constant case, when other 
data are combined, the allowed region lies around the same region as the ACDM case and 
also it is not so significantly larger compared to that for the ACDM model. Notice that 
the combination of CMB+SN again already gives severe constraint in the Qm^h plane. 
This result shows that even for a time-evloving equation of state, the SN data works well 
to break the degeneracy among and dark energy parameters for the parametrization 
Eq. (§. 

In Ref. [39], the so-called A parameter has been used to obtain a constraint from BAO, whereas Dy 
is used in this paper. However, as long as we also use the prior on ujm, these two constraints are almost 
equivalent. Also, the shift parameter R has been used to obtain a constraint from CMB but this gives a 
similar result to our method of combining 9a and uim- 
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0.1 0.2 0.3 0.4 0.2 0.3 0.4 0.5 



m ^ 

Figure 2: 2cr constraints are shown for the case with a constant equation of state for dark 
energy. We use the data from CMB alone (red sohd hne), CMB+BAO (green dashed 
hne), CMB+SN (blue dotted line) and all combined (purple shaded region). As for the 
treatment of SN data, we use those from Gold06 (left panel) and Davis07 (right panel). 
The value of Wq is marginalized over —3 < Wq < —1/3. 



We also performed the same analysis using the parametrization given in Eq. Q with 
= 0.5 and the results are shown in Fig. |4} As seen from the figure, the constraints on 
Qm and h are almost the same even if we assume this type of dark energy parametrization. 
Furthermore it should be noticed that here again the combination of CMB+SN gives a 
severe constraint. 

Although we have adopted just two types of parametrization of dark energy equation 
of state, as far as we consider some typical types of dark energy evolution, it seems that 
dark energy models do not affect much the determination of the Hubble constant using 
combined data of CMB, BAO and SN. As already mentioned and seen from Figs. [2| [3]and 
|4| CMB+SN gives sufficiently tight constraints in the flm^h plane. Meanwhile, the allowed 
region from CMB+BAO is rather broad and BAO does not seem to have much constraining 
power. (However, when we allow a non-flat universe, BAO becomes important, which we 



are going to discuss in Sec. 3.1.4 ) Thus we focus on the constraints from CMB and SN 
for a while and explain how the constraints can be obtained. 

Notice that since h and Qm are degenerate with respect to 9a, to break the degeneracy, 
we assume the prior on uJrn determined from the height of the CMB power spectrum. 
Thus the fact that h is well determined even if we allow some variations of dark energy 
models eventually indicates that the determination of Qm is not affected much by the 
assumption for dark energy. (Also notice that, as already mentioned, the SN data cannot 
determine the Hubble constant since the dependence on h in the luminosity distance is 
totally indistinguishable from the uncertainty in the absolute magnitude of SN.) Remind 
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Figure 3: Same as Fig. [2] except we assume the dark energy equation of state as Eq. ([T]) and 
marginalize over the values of Wq and Wi in the ranges — 3 < < — 1/3 and — 3 < Wi < 3 
with the prior Eq. 




Figure 4: Same as Fig. |2] except we assume the dark energy equation of state as Eq. Q 
with 2;* = 0.5 and marginalize over the values of wo and wi in the ranges —3 <wq < —1/3 
and — 3 < t^i < —1/3. 
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Table 3: The central values and la errors, along with the minimum value of from 
CMB+BAO+SN for the cases with various dark energy models. Here the dark energy 
parameters are marginalized over. The cases where the equation of state is always larger 
than —1 are also presented. A flat universe is assumed. The best-fit values of dark energy 
parameters are also shown. 
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that the SN data can give information of the background evolution from the present up 
to which represents the redshift of the furthest SN currently available (in observations 
used in the anlysis, zsn ~ 1-8). In this period, the background evolution is determined 
by dark energy and matter. Since the distance measure such as the luminosity distance 
involves some integration with respect to redshift, it is well-known that there is a severe 
degeneracy among Q^n and dark energy parameters. However we can break the degeneracy 
by combining the constraint from 6a which can also probe the background evolution earlier 
than zsN. In the epoch earlier than zsn? we can approximate that the universe is dominated 
by matter component. Thus, knowing the distance between the present epoch and ^sn by 
the SN data, the matter density is determined accurately from the distance measurement 
between zstss and z^-cc- The point is that the distance between z = and zsn, where both 
of the matter and dark energy contribute, can be inferred by the SN data independently 
from dark energy model. This is why the combination of SN and CMB can determine 
Qrn in a dark energy model-independent manner. We will provide a more quantitative 



demonstration of this point in Appendix B.l 



As we have just argued, the Hubble constant derived from CMB, BAO and SN is not 
affected much by the assumption for a dark energy model. Notice that, however, when we 
impose a prior on the Hubble constant by using a value for instance of the cosmic distance 
ladder measurement mentioned in Sec. 2.4[ the constraints on dark energy parameters 
can be affected by the prior on Hq. Naturally, its effect depends on the accuracy of the 
determination of the Hubble constant. We will return this issue in Sec. 13.21 

We summarize our results on the constraints on h and Qm for time-evolving dark 
energy equations of state in Table |3] along with the best-fit values of the dark energy 
parameters. As discussed above, even if we assume a different dark energy parametrization, 
the constraints on h and Qrn are almost unchanged. Furthermore, even if we restrict 
ourselves to the dark energy model with its equation of state being larger than —1 in 
the course of the history of the universe, the results are not affected by this assumption 
as seen from Table [sj It should be mentioned that two SN data sets (i.e., Gold06 and 
Davis07) give somewhat different values of the Hubble constant as discussed above. The 
difference due to the SN data set is more significant than the one caused by the dark 
energy parametrization. 



3.1.4 Case with a non-fiat universe 

Here we discuss the constraint on the Hubble constant without assuming a flat universe. 
When one considers the constraint on the Hubble constant from cosmological observations, 
a fiat universe is often assumed. However, to make our analysis general, here we allow a 
non-flat universe to obtain the constraint on the Hubble constant. The constraint on the 
curvature of the universe itself with some priors on the Hubble constant will be discussed 
in Sec. 13.31 

In Fig. |5| we show the constraints in the ^Im^h plane for the case with a cosmological 
constant while marginalizing over in the range —0.5 < fi^ < 0.5. We show the con- 
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Figure 5: 2cr constraints in the Vtm-h plane without assuming a flat universe. The cur- 
vature is marginalized over in the range —0.5 < fifc < 0.5. A cosmological constant is 
assumed for dark energy. Constraints from the data from CMB alone (red solid line), 
CMB+BAO (green dashed line), CMB+SN (blue dotted line) and all combined (purple 
shaded region) are shown. As for the treatment of the SN data, we use those from Gold06 
(left panel) and Davis07 (right panel). 
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Figure 6: Same as Fig. [s] (i.e. the curvature is marginalized over) except that a constant 
equation of state is assumed and marginalized over in the range —3 < wq < —1/3. 
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Table 4: The central values and la errors, along with the minimum value of from 
CMB±BAO±SN for the cases with various dark energy models. Here the dark energy 
parameters and Qk are marginalized over. The best-fit dark energy parameters and Qk are 
also shown. 



straints from CMB alone (red solid line), CMB±BAO (green dashed line), CMB±SN (blue 
dotted line) and all combined (purple shaded region) separately. The case with marginal- 
izing over a constant equation of state in addition is presented in Fig. [6] When we allow a 
non-flat universe, there is a strong degeneracy in the distance measures as in the case we 
assume a dynamical dark energy in a fiat universe which can be seen in Figs. [2| |3] and |4} 
However, when all the data are combined, there is little difference between the constraints 
for the cases with a non-fiat universe being allowed and a cosmological constant with a fiat 
universe. This is because the combination of all three cosmological data sets work well to 
remove the degeneracy even if we consider the possibilities of a non-fiat universe. On the 
other hand, when one focuses on the combination of two data sets such as CMB±BAO 
and CMB±SN, we can see some differences between a fiat case and a non-fiat case. Recall 
that when one assumes a fiat universe, the combination of CMB and SN already gives a 
relatively severe constraint even for the case with time-varying equations of state but that 
of CMB±BAO is not so severe. For the case with a non-fiat universe, in contrast to the 
fiat universe case, the combination of CMB±SN does not give a severe constraint. When 
one allows the nonzero curvature of the universe, the degeneracy in CMB (here this means 
9a) becomes worse: it involves flm, ^k and dark energy parameters. Thus just adding SN 
data is not enough to break this degeneracy. BAO is necessary to remove such degeneracy. 
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In Appendix [B| we discuss this point in a quantitative way. 

We have seen that the distance measures suffer from the degeneracies among param- 
eters describing dark energy equation of state and the curvature of the universe. Due to 
those degeneracies, if we use CMB alone, CMB+BAO and CMB+SN, we get different 
constraints as clearly seen in Figs. |5]and[6] It seems to be necessary to combine all three 
observations to remove the degeneracies. Once that is done, the constraints are almost the 
same as the case with the flat ACDM model. In Table |4| the central values and la errors of 
h and Qm are summarized for the cases with a cosmological constant, a constant equation 
of state, dark energy parametrizations given in Eqs. ([T| and ^ when a non-flat universe 
is allowed in the analysis. The table shows that the constraint on h is almost unchanged 
even if we allow a non-flat universe under some types of dark energy parametrization. 



3.1.5 Summary of constraints on the Hubble constant 

In Fig. [7} we summarize some of the results for the constraints on h and presented 
in this section. See also Tables |3] and |4] for more detailed constraints on h and flm and 
the best flt values of dark energy parameters and/or the curvature. As discussed above, 
even if we assume different types of dark energy parametrization, there is no considerable 
change in the constraints on h and flm, although the allowed region for the case with a 
cosmological constant and a flat universe is slightly smaller compared to the other cases. It 



is interesting to compare this result with the distance ladder measurements of h (Sec. 2.4). 
In particular, when all the data are combined, the central values of h cannot be as low as 
the Sandage's central value 0.62 and cannot be as high as the Macri's central value 0.74 
even if we relax the assumptions of a cosmological constant as dark energy and/or the 
flatness of the universe. Of course, since the current measurements of h by the distance 
ladder have somewhat large systematic errors (see Table [T]), this must be taken only 
as a quick comparison. The most seemingly problematic case arises when we compare 
CMB+BAO+Davis07 and Sandage's h, but their 1 a errors overlap. Our analysis shows 
that the allowed region of the Hubble constant from CMB, BAO and SN can be different 
from the WMAP flat ACDM value of h = 0.73 ± 0.03 by more than la depending on 
the assumption for dark energy and the cosmic curvature, but h < 0.59 or h > 0.76 are 
not allowed at 2a level which can be read off from Tables [3] and HI This conclusion is 
obtained for rather limited types of dark energy parametrization but the interpretation 
presented in the appendix leads us to speculate that this holds true for any dark energy 
parametrization. 

It should also be noticed that the choice of the SN data set makes larger difference than 
that of the assumptions on dark energy and the curvature of the universe. The Gold06 
data tends to give lower h and higher Qrn than the Davis07 data. Thus, there is much 
room for improvement in SN data as well as the distance ladder measurement of h. We 
need more precise measurement from both flelds of observation to tell whether there is a 
discrepancy or not. 

The errors on the Hubble constant (and on Q^) are slightly larger for the cases with 
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Figure 7: Contours Icr and 2a allowed regions in the Qm^h plane for the cases with: 
(a) a cosmological constant (same as Fig. [l] (d)), (b) a constant equation of state wx, 
marginalized over wx, (c) the parametrization defined by Eq. ([T]), marginalized over wq 
and Wi, and (d) the parametrization defined by Eq. Q with 2;* = 0.5, marginalized over 
Wo and Wi. Constraints using the SN data from GoldOG (red solid line) and Davis07 
(blue dashed line) are shown separately. In the panels (e)-(h), we marginalize over in 
addition to the panels (a)-(d) respectively. 
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a non-flat universe than those with a flat universe which can be seen from Tables [3] and 
m For a flat case, la errors on h are about 0.02 and, for a non-flat case, it is around 0.03. 
Therefore, when we use the current observations of CMB, BAO and SN, we can expect that 
this level of accuracy is needed for the distance ladder measurement of the Hubble constant 
to give a meaningful external prior on h in constraining the dark energy parameters and 
the curvature of the universe. Assuming that such sensitivity on h is achieved in the 
future, we investigate its effects on the determination of dark energy equation of state and 



the curvature of the universe in the following sections 3.2 and 3.3 



3.2 Constraint on dark energy equations of state 

Here we discuss the constraints on parameters of dark energy which describe the time 
dependence of its equation of state. Our special emphasis is on investigating how they are 
affected by several external priors on the Hubble constant and/or by the assumption of 
the flat universe. In most analyses on this issue so far, a flat universe is usually assumed. 
Once one invokes the inflationary paradigm, the flatness assumption seems to be natural. 
However, the flatness itself should be tested including the uncertainties of dark energy. In 
this respect, when we investigate the equation of state for dark energy, we also remove 
the assumption of a flat universe to obtain a more conservative constraint. In particular, 
in the following, we discuss implications of the results taking a cosmological constant as 
a reference point. Since a cosmological constant is the simplest and most conventional 
model for dark energy (and it can fit the observations satisfactorily as shown below), it is 
useful for illustrating the roles of the priors. 

First we assume the dark energy parametrization of Eq. ([T]). In Fig. |8| the constraints 
in the wq-wi plane are shown for several cases. In the panels (a)-(c), we marginalized over 
the values of Qm and h. In the panels (d)-(f), we repeat the same analysis except that we 
allow a non-flat universe and marginalize over flk- The best fit values of wo and wi and 
the minimum are summarized in Tables [s] and [g] The best fit values of h and (and 
ilk for the case of a non-flat universe in Table [6]) are also shown. In the panels (b), (c), 
(e) and (f), we assume some priors on the Hubble constant to see how it can affect the 
determination of dark energy parameters. We impose Gaussian priors as h = 0.72 ± 0.02 
(panels (b) and (e)) and h = 0.62 ±0.02 (panels (c) and (f)) where we use the Freedman's 
and Sandage's (see Table [T| as the central values and hypothetical errors of 0.02. 

The hypothetical errors adopted here are motivated from the results we obtained in 
Sec. 3.1 As summarized in Sec. [3.1.5 the uncertainties of h from the cosmological ob- 



servations are 0.02 ~ 0.03. Thus we can expect that this level of accuracy on the Hubble 
constant prior is required to have an influence on constraining dark energy sector. In 
fact, we have also made the analysis adopting a Gaussian prior with the current errors 
as shown in Table [T] but we cannot see a noticeable difference compared with the case 



of no prior on the Hubble constant. It should be noted that, as mentioned in Sec. 2.4 
this level of improvement on the Hubble constant determination is not unimaginable in 
near future [34]. Note that we can expect that other observations (CMB, BAO and SN) 
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Figure 8: la and 2a constraints from CMB+BAO±SN in the Wq-Wi plane marginalizing 
over and h are shown for the cases with (a) no prior on the Hubble constant, (b) 
assuming a Gaussian prior on the Hubble constant h = 0.72 ±0.02 and (c) h = 0.62 ±0.02. 
In the panels (d)-(f), we allow a non-flat universe and marginalize over Qk in addition to 
ilm and h. The black dashed lines show the boundary of the prior Eq. (|3]). The constraints 
using the SN data sets from Gold06 (red solid line) and Davis07 (blue dashed line) are 
shown separately. 
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Table 5: The best fit values for fi^, h, wo and wi for the analysis presented in the panel 
(a)-(c) of Fig. M The minimum values of are also shown. 
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Table 6: The best fit values for fi^; h, fi^, wo and Wi for the analysis presented in the 
panel (d)-(f) of Fig. M The minimum values of are also shown. 



would be more precise at the time when this level of accuracy in the Hubble constant 
determination is realized. Thus, our combining the hypothetical Hubble priors and the 
current observational data is not a forecast of future status in a strict sense but should 
rather be regarded as an illustration of how the Hubble external priors affect conclusions 
on properties of dark energy and the curvature of the universe. 

Let us start with looking at the constraints from the present cosmological observations 
assuming no Hubble prior. When we assume a fiat universe, although the best fit value 
is slightly away from a cosmological constant (especially in the positive direction of Wi. 
See Table |5]), a cosmological constant is within the 2a allowed regions for both SN data 
sets as shown in Fig. [s] (a). This situation holds when we allow a non-fiat universe and 
marginalize over ilk as shown in Fig. [s] (d). We can understand this by noticing that the 
region closely around a flat universe is favored in this parametrization with no Hubble 



prior as seen in Tabled or in the analysis presented in the following Sec. 3.3 



We now turn to the constraints when we assume the external priors on the Hubble 
constant. The constraints in the Wo-Wi plane for this case are shown in Fig. [s] (b), (c), (e) 
and (f). Probably the most interesting feature is seen in the panel (c): in a fiat universe, 
if the Sandage's value is confirmed at this level, a cosmological constant {wq = —I and 
Wi = in the figure) would be rejected at nearly 2a level in this parametrization. However, 
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as shown in the panel (f), if we allow a non-flat universe and marginalize over Qk, a 
cosmological constant is well within the allowed regions. This is an example which shows 
the importance of the priors on the Hubble constant and the curvature of the universe 
when we probe the nature of dark energy. 

Next we show the constraints on the dark energy parameters for the parametrization 
Eq. Q in Figs. |9] and 10 Although this parametrization includes three parameters, we 
present our result in the Wo^Wi plane fixing to several values: 0.1, 0.2, 0.5, 1.0 and 2.0. 
Namely we regard 2^, as labeling the model expressed as Eq. ^ which has two parameters 
Wq and Wi, and do not marginalize over z^,. In Fig. |9l a fiat universe is assumed and in 



Fig. 10, we allow a non-fiat universe and marginalize over the curvature. The best fit 
parameter values are shown respectively in Tables [7] and [8j We also impose some priors 
on the Hubble constant as is done for the analysis of the dark energy parametrization of 
Eq. 0. 



For the case with no Hubble prior, as the panels in the left column in Figs. M and 10 



indicate, a cosmological constant {wq = — 1 and wi = —1) is in the allowed regions with 
or without the assumption of a fiat universe. This is because, as is the case with the 
parametrization Eq. ([T]) discussed above, the region around the fiat universe is more or 
less favored in this parametrization with no Hubble prior as seen in Tables [7] and |8} 

However, if we impose some external priors on the Hubble constant, the constraints 
in the wo-wi plane receive some interesting effects of the external priors on the Hubble 
constant and the assumption for the curvature of the universe. We can see them from the 



panels in the middle and right columns in Figs. M and 10 



For example, when we impose h = 0.62 ± 0.02 for the z^: = 1.0 model assuming a 
fiat universe, a cosmological constant is not within the 2a allowed region for both SN 
data sets (see Fig. IqI). However, if we allow a non-fiat universe and marginalize over 



ilk, a cosmological constant is allowed with 2a level (see Fig. 10). This is similar to the 
tendency which is seen in the case of the parametrization Eq. (jl ) discussed above. Another 
interesting point can be observed in this parametrization, which is actually an opposite 
tendency to the one just mentioned. When we impose h = 0.72 ± 0.02 for the z^, = 0.5 
model assuming a flat universe, a cosmological constant is well within the allowed regions 
for both SN data sets (see Fig. [9]). However, when we use the GoldOG data set, if we allow a 
non-flat universe and marginalize over Qk, a cosmological constant is out of the 2a allowed 
region (see Fig. 10). This is considered to be an infrequent case in which relatively large 
Qk gives a significantly better fit than a flat universe (Xmin decreases by 4.3 in this case) 
with the equation of state well away from a cosmological constant, but we should bear in 
mind that such could happen. 

In summary, we have seen that the prior on the Hubble constant, the way to parametrize 
the equation of state and the assumption for the curvature of the universe are all very 
important to probe the nature of dark energy. These would be true for data from future 
experiments of CMB, BAO and SN. In particular, since there is a well-known severe de- 
generacy among wx and Qx, the independent determination of the Hubble constant will 
still affect the constraint on dark energy parameters by pinning down the value of Qm 
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Figure 9: la and 2a constraints from CMB+BAO+SN in the wq-wi plane marginalizing 
over h and Qm are shown for several values of z*. For each z^, we show the cases with no 
prior on the Hubble constant and Gaussian priors h = 0.72 ± 0.02 and h = 0.62 =b 0.02. A 
fiat universe is assumed. The constraints using the SN data sets from GoldOG (red solid 
line) and DavisOT (blue dashed line) are shown separately. 
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Figure 10: la and 2a constraints from CMB+BAO+SN in the Wq-Wi plane marginalizing 
over h, Q^n and flk are shown for several values of z^, for the cases with no Hubble prior 
and assuming Gaussian priors h = 0.72 ± 0.02 and h = 0.62 ± 0.02. 
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CMB+BAO+Gold06 


A mm 




h 


Wo 


Wi 


= 0.1 (No prior) 


157.9 


0.279 


0.678 


-2.52 


-0.71 


Prior h = 0.72 ± 0.02 


159.8 


0.257 


0.702 


-2.83 


-0.73 


Prior h = 0.62 ± 0.02 


161.2 


0.307 


0.648 


-2.11 


-0.68 


= 0.2 (No prior) 


157.5 


0.289 


0.665 


-1.63 


-0.68 


Prior h = 0.72 ± 0.02 


161.0 


0.262 


0.695 


-1.63 


-0.75 


Prior h = 0.62 ± 0.02 


159.6 


0.310 


0.643 


-1.65 


-0.62 


= 0.5 (No prior) 


157.5 


0.298 


0.654 


-1.21 


-0.61 


Prior h = 0.72 ± 0.02 


162.2 


0.265 


0.692 


-1.15 


-0.75 


Prior h = 0.62 ± 0.02 


158.6 


0.315 


0.637 


-1.24 


-0.54 


z^: = 1.0 (No prior) 


158.1 


0.304 


0.649 


-1.08 


-0.52 


Prior h = 0.72 ± 0.02 


163.1 


0.266 


0.692 


-1.01 


-0.77 


Prior h = 0.62 ± 0.02 


158.8 


0.319 


0.634 


-1.11 


-0.44 


z^ = 2.0 (No prior) 


158.6 


0.306 


0.649 


-1.02 


-0.38 


Prior h = 0.72 ± 0.02 


163.4 


0.266 


0.693 


-0.96 


-0.77 


Prior h = 0.62 ± 0.02 


159.3 


0.319 


0.635 


-1.01 


-0.33 


CMB+BAO+Davis07 


A mm 




h 


Wo 


Wi 


2* = 0.1 (No prior) 


195.5 


0.258 


0.707 


-2.00 


-0.85 


Prior h = 0.72 ± 0.02 


195.7 


0.251 


0.714 


-2.07 


-0.85 


Prior h = 0.62 ± 0.02 


204.1 


0.297 


0.665 


-1.76 


-0.80 


2* = 0.2 (No prior) 


195.9 


0.263 


0.700 


-1.40 


-0.84 


Prior h = 0.72 ± 0.02 


196.4 


0.253 


0.710 


-1.39 


-0.86 


Prior h = 0.62 ± 0.02 


203.2 


0.297 


0.662 


-1.54 


-0.72 


z^ = 0.5 (No prior) 


195.9 


0.267 


0.692 


-1.18 


-0.77 


Prior h = 0.72 ± 0.02 


196.7 


0.255 


0.708 


-1.14 


-0.84 


Prior h = 0.62 ± 0.02 


201.2 


0.298 


0.653 


-1.31 


-0.57 


2* = 1.0 (No prior) 


195.5 


0.273 


0.684 


-1.15 


-0.61 


Prior h = 0.72 ± 0.02 


196.8 


0.256 


0.706 


-1.09 


-0.77 


Prior h = 0.62 ± 0.02 


199.0 


0.303 


0.646 


-1.24 


-0.38 


= 2.0 (No prior) 


195.4 


0.277 


0.678 


-1.12 


-0.38 


Prior h = 0.72 ± 0.02 


196.9 


0.257 


0.706 


-1.06 


-0.67 


Prior h = 0.62 ± 0.02 


199.3 


0.301 


0.653 


-1.08 


-0.33 



Table 7: The best fit values for flm, h, wo and wi for the analysis presented in Fig. |9] in 
which a fiat universe is assumed. The minimum values of are also shown. 
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CMB+BAO+Gold06 


Amm 




h 




Wo 


Wi 


= 0.1 (No prior) 


157.8 


0.273 


0.684 


0.003 


-2.66 


-0.68 


Prior h = 0.72 ± 0.02 


158.4 


0.251 


0.714 


0.013 


-3.20 


-0.60 


Prior h = 0.62 ± 0.02 


159.7 


0.318 


0.636 


-0.013 


-1.58 


-0.85 


= 0.2 (No prior) 


157.3 


0.278 


0.678 


0.009 


-1.80 


-0.59 


Prior h = 0.72 ± 0.02 


158.3 


0.253 


0.711 


0.025 


-2.06 


-0.49 


Prior h = 0.62 ± 0.02 


159.2 


0.316 


0.637 


-0.010 


-1.44 


-0.75 


z^, = 0.5 (No prior) 


156.6 


0.279 


0.676 


0.032 


-1.39 


-0.33 


Prior h = 0.72 ± 0.02 


157.9 


0.255 


0.709 


0.042 


-1.38 


-0.33 


Prior h = 0.62 ± 0.02 


158.6 


0.314 


0.638 


0.002 


-1.25 


-0.51 


z^: = 1.0 (No prior) 


157.5 


0.294 


0.661 


0.018 


-1.12 


-0.33 


Prior h = 0.72 ± 0.02 


160.0 


0.260 


0.704 


0.032 


-1.09 


-0.33 


Prior h = 0.62 ± 0.02 


158.7 


0.318 


0.635 


0.009 


-1.14 


-0.33 


z^ = 2.0 (No prior) 


158.5 


0.302 


0.653 


0.005 


-1.01 


-0.33 


Prior h = 0.72 ± 0.02 


161.7 


0.263 


0.702 


0.020 


-0.98 


-0.33 


Prior h = 0.62 ± 0.02 


159.3 


0.321 


0.633 


-0.003 


-1.03 


-0.33 


CMB+BAO+Davis07 


A mill 




h 




Wo 


Wi 


2* = 0.1 (No prior) 


195.1 


0.271 


0.686 


-0.008 


-1.70 


-0.93 


Prior h = 0.72 ± 0.02 


195.7 


0.251 


0.713 


-0.001 


-2.04 


-0.86 


Prior h = 0.62 ± 0.02 


197.5 


0.312 


0.639 


-0.023 


-0.94 


-1.09 


2* = 0.2 (No prior) 


195.4 


0.278 


0.678 


-0.010 


-1.20 


-0.97 


Prior h = 0.72 ± 0.02 


196.4 


0.253 


0.711 


0.000 


-1.38 


-0.86 


Prior h = 0.62 ± 0.02 


197.4 


0.313 


0.638 


-0.023 


-0.94 


-1.14 


2* = 0.5 (No prior) 


195.5 


0.281 


0.675 


-0.011 


-1.07 


-0.98 


Prior h = 0.72 ± 0.02 


196.7 


0.254 


0.710 


0.002 


-1.16 


-0.79 


Prior h = 0.62 ± 0.02 


197.3 


0.315 


0.637 


-0.026 


-0.93 


-1.29 


z^: = 1.0 (No prior) 


195.4 


0.279 


0.677 


-0.008 


-1.09 


-0.84 


Prior h = 0.72 ± 0.02 


196.5 


0.253 


0.710 


0.013 


-1.17 


-0.47 


Prior h = 0.62 ± 0.02 


197.3 


0.315 


0.637 


-0.026 


-0.96 


-1.52 


z^, = 2.0 (No prior) 


195.4 


0.278 


0.677 


-0.005 


-1.10 


-0.53 


Prior h = 0.72 ± 0.02 


196.6 


0.255 


0.709 


0.011 


-1.10 


-0.33 


Prior h = 0.62 ± 0.02 


197.3 


0.315 


0.637 


-0.027 


-0.96 


-2.08 



Table 8: The best fit values for Qm, h, flk, wo and wi for the analysis presented in Fig. [TO 
in which a non-flat universe is allowed. The minimum values of are also shown. 
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more precisely. The answer to the simplest question of "a cosmological constant or not" 
can be altered by changing some of these assumptions. 



3.3 Constraint on the curvature of the universe 

Here we discuss the constraints on the curvature of the universe assuming some priors on 
the Hubble constant including some types of dark energy model in addition to a cosmolog- 
ical constant and a constant equation of state using the parametrizations Eqs. ([T]) and (|4]). 
As already mentioned, a flat universe is assumed in most cosmological parameter estima- 
tion in the literature because the inflationary paradigm strongly suggests it. However, the 
paradigm should be tested through the test of the flatness and it should be done includ- 
ing the uncertainty in dark energy sector since its nature is not understood yet. In this 
respect, the investigation of the curvature of the universe in connection with dark energy 
models have been done in Refs. [4,39-55]. In Refs. [39-41], it was shown that an open uni- 
verse can be largely allowed for some particular dark energy parametrizations. However, 
in the previous works, the Hubble constant was treated as a nuisance parameter to be 
marginalized over. Here we analyze how the constraint on the curvature is correlated with 
the Hubble constant and how the priors on it affect the constraint. We also investigate 
whether the constraint varies according to the choice of dark energy parametrization. 



In Fig. 11 , the constraints in the Qk^h plane are shown for the cases with a cosmological 
constant, a constant equation of state, time-varying equations of state paramtrized as 
Eq. ([T]) and Eq. ^ with = 0.5. For each model, we impose no prior and two types of 



prior on the Hubble constant as in Sec. 3.2 



From the current data of CMB, BAO and SN without an external Hubble prior, we 
may conclude that the universe is constrained to be around fiat. In almost all of the cases, 
the curvature is limited as \flk\ < 0.05 at 2a level (only exception is that the case where 
the Gold06 set is used for the parametrization Eq. Q with z^: = 0.5. The 2a boundary 
extends to flk ~ 0.07). However, when we look at them more closely, we notice some 
dependence on dark energy parametrization. For the cases with a cosmological constant 
and a constant equation of state, the allowed region extends into a closed universe, which 
is a rather widely known result (e.g. Fig. 17 of Ref. [1]). In contrast, for the case with a 
time-varying equation of state Eq. Q with z^ = 0.5, the allowed region has a much wider 
area in an open universe as has been found in Ref. [40]. 



Fig. 11 also shows the known tendency that the lower Hubble constant is favored in a 
closed universe (a positive correlation between h and Qk)- Due to this correlation, when 
we assume the prior oi h = 0.62 ± 0.02 for the Hubble constant, the region with a closed 
universe occupies larger space in the allowed region. On the other hand, when the prior 
h = 0.72 ± 0.02 is assumed, wider region of an open universe is allowed compared to a 
closed universe. For example, a fiat universe is rejected at 2a level and a closed universe 
is favored for a cosmological constant with the prior h = 0.62 ± 0.02. Meanwhile, the 
prior h = 0.72 ± 0.02 enhances the preference of an open universe to 2a level for the 
parametrization Eq. ^ with z^ = 0.5. 
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Figure 11: la and 2a constraints from CMB±BAO±SN in the Qk^h plane marginalizing 
over flm and dark energy parameters. The dark energy parametrization adopted in the 
analysis are shown on the right side of the panels. We show the cases with no prior on the 
Hubble constant and Gaussian priors h = 0.72 ±0.02 and h = 0.62 ±0.02. The constraints 
using the SN data sets from Gold06 (red solid line) and Davis07 (blue dashed line) are 
shown separately. 
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-0.004 


0.706 


0.258 


-1.01 
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197.6 
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0.311 


-1.07 




Wq and Wi 


195.4 


-0.009 


0.676 


0.280 


-1.10 


0.37 


Prior h = 0.72 ± 0.02 


196.6 


0.006 


0.710 


0.254 


-1.15 


0.81 


Prior h = 0.62 ± 0.02 


197.3 


-0.026 


0.638 


0.314 


-0.93 


-0.93 


Wq and wi = 0.5) 


195.5 


-0.010 


0.676 


0.280 


-1.07 


-0.99 


Prior h = 0.72± 0.02 


196.7 


0.002 


0.710 


0.254 


-1.15 


-0.80 


Prior h = 0.62 ± 0.02 


197.3 


-0.025 


0.638 


0.315 


-0.94 


-1.28 



Table 9: Best fit values of and h for various assumptions on dark energy. Here we also 
show the best fit values for marginalized parameters such as Qrn,WQ{wQ) and Wi{wi). 



The final remark is that, in the previous works [40,41], it has been shown that when 
we parametrize dark energy equation of state as Eq. Q and take 2;* ~ 0.5, the region of an 



open universe as large as ilk ~ 0.2 is allowed. However, for the result presented in Fig. 11 
this is not the case even though we use the same type of parametrizaion. This is because 
we in this paper use a strong prior wx < on the equation of state given in Eq. ^ 



but in the previous works [40,41], a weak prior wx < has been adoptee As is already 
noticed in Refs. [40,41], an open universe tends to be preferred for dark energy whose 
equation of state approaches to zero at earlier time. Since we remove this possibility by 



"^"^ The method of this paper to obtain a constraint from cosmological data is also slightly different 
from that used in [39-41]. However, the difference in the method is irrelevant to the conclusion on the 
constraint on the curvature of the universe. 
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adopting a strong prior, we do not have a large allowed region with an open universe. In 
other words, our result here confirmed explicitly that the dark energy whose equation of 
state is close to zero at earlier time, which is sometimes called early dark energy, can help 
to allow an open universe at least as long as the background evolution is concerned. 

4 Conclusions and Discussion 

We studied the constraints on the Hubble constant from CMB, BAO and SN assuming 
several types of dark energy parametrization. Although the Hubble constant and dark 
energy are both important issues in cosmology today, these two subjects have not been 
investigated much simultaneously. First we investigated the constraints in the flm^h plane 
assuming several dark energy parametrizations. The constraints on the Hubble constant 
from the combination of CMB, BAO and SN observations obtained under the different 
dark energy models and priors are summarized in Table [3] when a flat universe is assumed 
and in Table |4] when the flatness assumption is dropped respectively. It is noticed that 
the constraints are not affected drastically by the dark energy model assumed and/or the 
assumption of the flatness of the universe. It is rather more affected by the choice of the 
the SN data set: the Gold06 set gives slightly lower value of h than the Davis07 set. We 
can conservatively conclude that Hq < 59 and Hq > 76 are highly unlikely from these 
cosmological observations: these parameter regions are not allowed at 2a level for any 
dark energy parametrization even if we do not restrict ourselves to a flat universe. Since 
the distance ladder estimations of Hq have somewhat large systematic errors at present, 
we are not at the stage of arguing any possible discrepancies among these measurements 
now. Nevertheless, it is worth mentioning in passing that the Sandage's central value 
Hq = 62 [33] and the Macri's value Hq = 74 [34] are fairly close to the limit we have 
obtained here. 

We have also investigated the constraints on some dark energy parameters assuming 
several priors on the Hubble constant. The constraints are derived with and without the 
assumption of the flatness of the universe. Using the present cosmological observations 
assuming no Hubble prior, we have found that a cosmological constant and a flat universe 
can fit all the data satisfactorily. When we impose a prior on the Hubble constant, we have 
adopted Gaussian priors h = 0.72 ±0.02 and h = 0.62 ±0.02. The central values are those 
of Freedman's [32] and Sandage's [33] but the error are taken to be a hypothetical value 
to give a meaningful effect on the dark energy parameter estimation. We have found that, 
even with some limited options discussed in this paper for the assumptions with respect 
to the prior on the Hubble constant, a parametrization of the dark energy equation of 
state, the curvature of the universe, they can have a great influence on determining the 
nature of dark energy. It should also be mentioned that the choice of the SN data set 
affects the allowed region. We demonstrate these points by taking a cosmological constant 
as a reference dark energy model because it is the simplest and most conventional model 
which can fit the observational data sets. For example, when we adopt the parametrization 
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Eq. ([T]) in a flat universe, the prior h = 0.62 ± 0.02 reject a cosmological constant at 2a 
level for the Gold06 data set, but the allowed region broadens to allow a cosmological 
constant if we do not adopt the flatness assumption or if we instead adopt the prior 
h = 0.72 ± 0.02. Moreover, if we adopt the parametrization Eq. ^ with z^, = 0.5 and 
the prior h = 0.72 ± 0.02, the Gold06 data set is in good agreement with a cosmological 
constant in a flat universe, but if we drop the flatness condition and marginalize over the 
curvature, a cosmological constant is disfavored at 2a level. These examples imply that 
our understanding of the nature of dark energy can be varied by the assumptions on the 
Hubble constant, a parametrization of dark energy equation of state and the curvature of 
the universe. 

Finally, we have investigated the constraints on the curvature of the universe assum- 
ing several types of dark energy and the priors on the Hubble constant. In contrast to 
the constraints on the Hubble constant, we see the result depends on the dark energy 
parametrization we adopt. For the cases with a cosmological constant and a constant 
equation of state, the allowed region occupies larger area in a closed universe, whereas for 
the case with a time-varying equation of state parametrized as Eq. Q with = 0.5, the 
allowed region extends into a region of an open universe. Since there is an obvious positive 
correlation between h and Qk, the prior h = 0.62 ± 0.02 exaggerates the preference for a 
closed universe and the prior h = 0.72 ± 0.02 for an open universe. This is what we see 



in Fig. 11 We also reconflrmed that the preference of an open universe is caused by the 
equation of state for dark energy which is close to at earlier time. This is because, we 
have found that flk as large as 0.2 is allowed in our previous papers Refs. [40,41] under 
the weak prior wx < 0, whereas Qk is found to be well below 0.1 in this paper under the 
stronger prior Wx < —1/3. 

Since dark energy is one of the most important problems in science today, a large 
number of works are focusing on dark energy itself. However, when one tries to probe the 
nature of dark energy, other cosmological parameters such as the Hubble constant, which 
was discussed in this paper, should necessarily be involved in various manners. In light 
of precise measurements of cosmology that we are having now, the works from this kind 
of viewpoint should be done to check our understanding of cosmology and may also give 
insights to probe the present state and the evolution of the universe. 
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Appendix 



A Effects of dark energy perturbation on CMB 

In the analysis for the constraints, we make use of the quantities which can be thor- 
oughly determined by the background evolution. Thus we only take into account the the 
modification to the background evolution by dark energy to obtain constraints on some 
parameters although dark energy component can fluctuate in general to affect the comic 
density fluctuation such as CMB anisotropies. As mentioned in the text, when dark en- 
ergy becomes dominant component at late time, the effect of fluctuation is not significant 
except on large scales and that from the modification to the background is enough to 
extract the cosmological constraints. However, when the equation of state for dark energy 
approaches to zero at earlier time, which means that the energy density of dark energy 
behaves px oc and can be comparable to that of matter, fluctuation of dark energy 
becomes important to affect the structure of acoustic peaks. In such a case, fluctuation of 
dark energy should be properly taken into account in the analysis. Thus, as for the CMB, 
the constraint from the distance measure such as the acoustic scale which is given by the 
information on the background evolution alone becomes invalid. This is one of the reason 
why we assume the prior of Eq. ^ to avoid such a case where dark energy has significant 
fraction at earlier time. In addition, when fluctuation of dark energy becomes important, 
other nature of dark energy such as the effective sound speed, which we denote here as 
(?xi can also modify density fluctuation. In this sense, the equation of state is not enough 
to consider the effect of dark energy. 

To see how the nature of dark energy can affect CMB, we show the CMB power spectra 



in Fig. 12 for several cases of the equation of state and the effective sound speed. Here 



we assume the parametrization of Eq. ([T]) for dark energy In the left panel, we fix the 
value of wq as wq = —1 and vary wi as wi = 1 (red solid line), 0.8 (green dashed line) and 
(blue dotted line). For other cosmological parameters, we assume the mean values of a 
power-law ACDM model from WMAP3 alone analysis as Um = 0.1277, cu;, = 0.02229, r = 
0.089, h = 0.732 and Ug = 0.958. As seen from the figure, when we compare the cases with 
wi = and 0.6, the power spectra is just shifted to smaller I but the structure of acoustic 
peaks is unchanged. This shift is caused by the the change of the angular diameter distance 
to last scattering surface due to the modification to the background evolution. However, 
when we take i^i = 1 in which energy density of dark energy is not negligible compared 
to that of matter at earlier time, the structure of acoustic peaks is significantly modified 
because fluctuation of dark energy can affect it in addition to the shift of the peaks in 
this case. Furthermore, this change depends on the nature of dark energy fluctuation. 



In the right panel of Fig. 12, we plot the cases with cj^ = 1 (red solid line), 0.1 (green 



dashed line) and (blue dotted line) for a parameter set {wo,Wi) = (—1, 1). Even though 

Other parametrization which makes dark energy density non- negligible at the epoch of recombination 
is adopted in Ref. [56,57] and constraints on the dark energy parameters are investigated. 
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Figure 12: (Left panel) CMB power spectra for a dark energy model with the parametriza- 
tion of Eq. ([T]). Here we fix the value of Wq = —1 and take several values for wi as Wi = 1 
(red solid line), 0.6 (green dashed line) and (blue dotted line). In this panel, we assume 
the effective sound speed as Cx = 1. (Right panel) CMB power spectra for the equation of 
state wq = — 1 + (1 — a) (i.e., we take {wo,wi) = (—1, 1) for the parametrization of Eq. ([T]). 
) The cases with c\ = 1 (red solid line), 0.1 (green dashed line) and (blue dotted line) 
are shown. 



the equation of state is the same (namely, the background evolution is the same), the 
CMB power spectra are different when one assumes different sound speed. In this kind 
of case, since the effective sound speed (and/or possibly another property of dark energy 
fluctuation) can affect the CMB, one should take into account the whole information of 
CMB power spectrum and constraints from the background evolution become invalid. 

We made another plot. Fig. 13, to see when neglecting the fluctuation of dark energy 
component is valid. In the left panel, we plot the height of the CMB first peak as a 
function of Wi for the parametrization Eq. ([T]) with fixing wq = —\. We can see that it is 
unchanged for Wi ^ 0.7, which indicates that the fluctuation of dark energy does not affect 
the structure of the acoustic peaks. In the right panel, as a function of wi, we plot the ratio 
of energy densities of dark energy and matter at the recombination epoch, z^ec = 1089. 
Notice that when wi ^ 0.7, energy density of dark energy is negligible compared to that 
of matter at recombination epoch. In other words, the effect of fluctuation of dark energy 
can be neglected when energy density of dark energy is small enough at earlier time. In 
such a case, since the information on the background evolution alone well captures the 
effect of dark energy, we can constrain dark energy parameters by only studying the shift 
of acoustic peaks. Thus, our method to obtain constraints from observations of CMB is 
justified when the equation of state is in a range where the energy density of dark energy 
is negligible at earlier time. The prior we take in this paper, Eq. ^ wx < —1/3, can 



satisfy this requirement. The final comment is that, as can be inferred from Fig. 13 



slightly looser prior like wx < —0.3 may make dark energy subdominant at the epoch of 
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Figure 13: (Left panel) The height of the first peak as a function of wi of the parametriza- 
tion Eq. ([T]) for the cases with = 1 (red sohd hne) and (green dashed hne). The value 
of wq is fixed as wq = —1. (Right panel) The ratio of the energy density of dark energy to 
that of matter at z = 1089 function of Wi. 



recombination and justify our analysis (and would not change our results much). However, 
since this condition depends on other cosmological parameters such as fix and Qm, we 
adopt the conservative prior of wx < —1/3. 



B How and ^ire determined from CMB, BAO 
and SN 

In Sec. |3.1[ we found that h is constrained rather tightly regardless of the dark energy 
model and/or the assumption of the flatness of the universe if all the CMB, BAO and SN 
data are combined. We also argued that this is equivalent to the determination of from 
the acoustic scales, 9a and Dy(0.35), and the luminosity distance cIl since Um = fimh"^ is 
given by the height of the CMB acoustic peak. Thus, in this appendix, we explain how 
these three types of cosmological data sets can determine Qm without referring to any 
particular dark energy model and without assuming the fiat universe. 



B.l What SN and CMB determine 

In this section, we focus on SN and CMB observations. We start with defining 



where zsn denotes the highest redshift to which SN data can probe. For the present data, 
we take z<^-^ = 1.8. Using this integral, the comoving angular distance to the last scattering 
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m 



Figure 14: Contours of rg{zrec) = 14.3 Gpc for several values of Qk (black solid lines) and 
/sN derived from observations of Gold06 (red dotted line) and Davis07 (blue dot-dashed 
line) . 



surface Eq. (|6| is given by 




dz' 



3.0 X lO^Mpc 



^snJ 



-1/21 



} • 



(26) 
(27) 



In the second line, we analytically performed the integration by approximating the uni- 
verse to be matter-dominated for 2;sn < z < Zrec and neglected the term proportional to 
(1 + z^cc)~^^'^ since ^sn ^^rec- Now, if JsN is known from SN observation, since CMB 
observation gives Um and rg{zrec), we obtain the degeneracy relation between Qm and 
Qk- This is drawn in Fig. 14 as a contour of rg{z-rec) = 14.3 Gpc (this is derived from 
Eq. (16) and rs{z^ec) = 149 Mpc, which is in turn from Eq. (10) using Um = 0.1277 and 
Ub = 0.02229) in the fim-/sN plane for several values of flk- Specifically, if we impose a 
flat universe prior, Qk = 0, Qm is determined uniquely. In Ref. [45], essentially the same 
argument is used to investigate the prospect of measuring Qk very precisely. 

The next issue is how tightly Jsn is constrained by the data. The SN data such as 
Gold06 or Davis07 used in this paper consist of the z-fi relation with an error on n where 



IjL is given by Eq. (22). The integral /sn is directly obtained from /^(^sn) if Hq and M are 
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known. However, the uncertainties in them cannot be resolved by SN data alone. Thus to 
cancel these uncertainties, we use /i at another redshift Zi which is close to zero, {zi <^ 1 is 
required not to make dependent on cosmology. Notice that, as we will adopt below, 
approximately we can use the relation ^^(-z) ^ z/Hq for small z.) To cancel the constants, 
we take a difference of ii{z) as 

M^SN)-M^i) = 51og^f^. (28) 



Since 



d^izsn) = l^^ S (v^/sn) , (29) 



and, (Il^Zi) ^ Zi/Hq for 2;i ^ 1, we obtain 

JsN = -^S-' ( 5ivfMio(M^sN)-M^O)/5^ . (30) 

\ 1 + ^SN J 

We derive /^(zsn) and fi{zi) by fitting the SN data to a power-law function fi{z) = az^. 
For the Gold06 set, 

ac = 44.288 ± 0.022, (31) 

Pg = (6.048 ± 0.033) x 10"^ (32) 

and for the Davis07 set, 

ttD = 44.278 ± 0.023, (33) 

Pd = (6.123 ± 0.030) x 10"^ (34) 

Note that it does not refer to a specific parametrization for dark energy. The x^'s are 
158.9 for the Gold06 data and 199.0 for the Davis07 data which are comparable to the 
values by the usual approach [41,58] of assuming some dark energy models, /sn has an 
uncertainty as regards the value of zi. However, a numerical experiment reveals that /sn 
barely depends on the value of zi for 0.02 ^ zi^ 0.05 {z ~ 0.02 is the lowest redshift for 
the SN data we here use). Thus, we take Zi = 0.05 hereafter. Then in a fiat universe, 

JSN,G = 1.10, (35) 
/sN,D = 1.15, (36) 

respectively for Gold06 and Davis07. Notice that they do not depend much on the cur- 
vature. In fact, they change less than 1% for \Qf^\ < 0.2. These values are plotted as 
horizontal lines in Fig. [I4j They cross the contour which satisfies the CMB observation 
for a flat universe (denoted by the black thick solid line) at around flm = 0.25. This 
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is consistent with the full analysis result as presented in Sec. 3A_ Moreover, we can see 
from the figure that since /sn,g < -^sn,d (which in turn comes from the difference in the 
power-law slope pc < Pd), Gold06 favors larger flm than Davis07 does. This is another 
point which we have noted in Sec. |3.1 We believe that these explanations are independent 
of the dark energy model we adopt because we use the SN data without referring to dark 
energy. 

However, Fig. 14 also shows that if we abandon the flatness assumption and shift flk 
from zero, different value of Qm is favored. This is why we cannot constrain Qrn much 
from CMB+SN when we marginalize over Q^. 



B.2 What BAO adds 



The BAO data gives -Dy (-Zrao) defined by Eq. (19) where ^bao = 0.35. This is written as 

n / \3 /■ \2 rj-l zbao 

Dv{zbao) = reizBAo) Hq 



^/Qra{l + ZbAoY + ^fc(l + ZbAO^ + {1 - - ^^fc)/(^BAo) ' 

(37) 

where f{z) is a function which expresses the evolution of the dark energy density. Here, 
the dark energy dependence enters in two places, re^ZBAo) and /(zbao)- Since zbao is 
relatively small, we may neglect the dark energy evolution (namely, we approximate dark 
energy as a cosmological constant) to approximate as f{zBAo) ~ 1- For '"^(^bao), since 
this is written as 



?(^bao) 



-^bao 



(3^ 



where /rao is the integral similar to Eq. (25 ) with zsn replaced by zbao, it can be calculated 



without referring to a dark energy model provided that we combine with the SN data. This 
is because we can infer /rao by using the fit to the SN data in the same manner to derive 



JsN in Appendix B.l Similarly to Jsn, -^bao depends slightly on the SN data (but not on 
Zi or the curvature). Gold06 and DavisO? give respectively 



-^BAO,G 
-^BAO.D 



0.31, 
0.32. 



(39) 
(40) 



Thus, here and in Appendix |B.l 
of CMB and BAO observables (ro(z. 



Now, if we fix Um using the CMB value, we have a relation between and determined 
from SN and BAO data for a measured -Dy(0.35). 

we have replaced the dark energy dependent part 
■) and Dy(0.35)) by the empirical values (Jsn and 
Ibao) inferred from the SN data. Namely, we can now draw contours of re{zrec) and 
Dy(0.35) in the fi^-fifc plane without referring to any dark energy model, 
we draw the contours rQ^Zj-ec 
Dy(0.35) 



In Fig 15 



14.3 Gpc (see Appendix B.l) by the black solid line and 
1402 Mpc (see Sec. 2.2) by the black dotted line. We have fixed cUrn to the 
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CMB value of 0.1277. Since the two SN data sets give slightly different values of /sn and 
Jbao as mentioned above, we show the results in two panels separately (the Gold06 data 
are used in the left panel and the Davis07 data in the right panel). The rQ^Zj-ec) contours 
run in somewhat diagonal direction, showing the degeneracy between and Qk for the 



CMB+SN combination as mentioned in Appendix B.l In contrast, the Dy(0.35) contours 



run almost horizontally, showing that the BAO+SN combination is insensitive to Qk and 
can determine flm regardless of the assumption on the curvature of the universe. This 
is reasonable because BAO measures the distance to relatively low redshift {z = 0.35). 
Ref. [25] has derived a linearized relation from BAO measurement: flm = 0.273 + 0.137^^ 
(Eq. (6) in Ref. [25]), which is quite similar to our slope for the Dv{0.35) contours. But also 
note that this relation has been derived only for a cosmological constant. If we consider 
the case with a constant equation of state or a time- varying equation of state, Qm cannot 
be determined from BAO alone. Our point is that when BAO and SN are combined, Qm 
can be determined regardless of the assumption on the curvature of the universe and dark 
energy model. 

In order to check the validity of these approximations, we plot the allowed regions by 
the full analysis such as explained in Sec. [2] The re(zrec) contours are compared with the 
analysis using constructed from 9a, tJm and the SN data set. The Dy(0.35) contours 
are compared with the one using from Dy(0.35), Um and SN. For completeness, we 
also performed the all-combined analysis using 6a, -Dy(0.35), Um and SN, which can be 



compared with the analyses of Qm and Qk done in Sec. 3.1| and Sec. 3^ The x^'s are 
minimized over h and dark energy parameters. We analyzed several dark energy models 
parametrized as Eqs. ([T]) and Q but the results do not show much difference. In FigJTsj we 
show 2a allowed regions for the cases with Eq. ([T| (blue dot-dashed lines) and Eq. (|4j) with 
z^: = 0.5 (red dashed lines) as examples. Note that the allowed regions for two models look 
alike. We can also see that the regions basically extend in the directions of the approximate 
re{z^cc) and DviO.35) contours. Therefore, it is considered to be appropriate to drop the 
dark energy dependence by using the SN data as we have done here. 
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Figure 15: Contours of rg{z,.cc) = 14308.5 Mpc (Eq. (27)) are drawn in black solid lines and 
those of (0.35) = 1402 Mpc are drawn in black dotted lines. Here, re^z^ec) and Dv{0.35) 
are approximate expressions in which dependence on dark energy model is dropped by 
using empirical values from SN data as discussed in Appendix |B.1| and |B.2[ Also shown 



are 2a allowed regions derived by the full analysis assuming some dark energy models 
and marginalizing over two dark energy parameters and h. The models are: (i) the 
parametrization Eq. ([T| and (ii) the parametrization Eq. Q with z^: = 0.5. The bands 
running diagonally are from SN, 6a and Um (blue dot-long-dashed line for Model (i) and 
red long-dashed line for Model (ii)), and the horizontal bands are from SN, Dv{0-35) and 
Um (blue dot-short-dashed line for Model (i) and red short-dashed line for Model (ii)). 
The shaded regions are from all data combined, SN, 6a, Dv{0.35) and (the regions 
are bounded by blue dot-dashed hne for Model (i) and red dashed line for Model (ii)). 
We can see that the contours of the approximate expressions well indicate the degeneracy 
directions of the corresponding data combinations. We use the Gold06 SN data set in the 
left panel and the Davis07 data in the right panel. 
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